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It has been recently proposed to use the operator product expansion to evaluate the expectation values of renor- 
malized operators without the need of a direct computation of the relevant renormalization constants. We test the 
viability of this idea in the two-dimensional non-linear <r-model discussing the non-perturbative renormalization 
of the energy-momentum tensor. 



The Operator Product Expansion (OPE) con- 
sists in rewriting, for y — ► x, 



A(x)B(y) 
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x - y)C{y) 



(1) 



Here A, B, and C are operators; Wq (x — y) 
are c- number functions (Wilson coefficients), sin- 
gular when | a; — y\ — > 0, that can be computed, 
for example, in perturbation theory. Dimensional 
analysis tells us that the leading contribution for 
y — > x is due to the operators C's in Eq. ([!]) 
that have the lowest dimension. Notice that the 
OPE is an operator relation and therefore, for any 
matrix element (ip\A(x)B(y)\4>), the coefficients 
Wq(x — y) are independent of the states 
and \cf>). This also means that if one wants to use 
the OPE, one should go through all the subtleties 
in the definition of the operators in quantum field 
theory (regularization, renormalization, . . . ). 

The authors of reference have proposed to 
consider, for lattice applications, the OPE in the 
particular case in which A and B are the com- 
ponents of a conserved current J M . This case is 
particularly simple because the components of J M 
do not need to be renormalized. This means that 

[Mx)j v {o)] MS (m) = [Mx)j v (o)] L (Vo) . ( 2 ) 

where the currents are expressed in terms of the 
bare fields and the coupling constant on the lat- 
tice is related to the coupling constant g in the 
continuum at the scale fi by a renormalization 
factor 



g L {l/a) = Z% (g(fi) ) fia) g(n) . 



(3) 



Therefore, the lattice computation of 
{tp\ Jfj, (x) J v (0)| <j>) provides the same expectation 
value as in the continuum, say in the MS-scheme, 
at the scale One can then use Eq. (|l|) to es- 
timate the matrix elements of the operators C's 
that dominate the r.h.s. of (Q) (in the continuum 
at the scale fi), bypassing all the troubles of lat- 
tice perturbation theory (large deviations from 1 
of the if's, large corrections, improvements, . . . ). 
The recipe is the following. One evaluates on 
the lattice the matrix elements (tf}\J fi (x)J v (0)\(j)) 
in the range 1 < \x\ < £ (where £ is as usual 
the correlation length) and then uses the OPE 
in the continuum (at the scale fj,) to fit the 
(now) numbers (ip\C(Q)\<f>). One needs the Wil- 
son coefficients Wq V {x). The idea is to use their 
expression computed in perturbation theory in 
the continuum. 

Assuming C(0) to be multiplicatively renor- 
malized, one can then use 

^\C(Q)\4>)=Zc{gL{l/a)^a)^\C L m4>) ■ (4) 

In this way, by measuring (ip\C L (0)\(f)) on the lat- 
tice, we obtain a non-perturbative determination 

oiz c . 

The discussion in Q ends with a question: is 
it feasible? That is, does it exist, in practice, a 
numerically accessible window for |jc| where, for 
a truncation of the sum over C, the OPE applies 
accurately enough to extract the expectation val- 
ues of the C's? 

To perform a check we need a case in which we 
already know the matrix elements (i/>|C(O)|0), so 
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that we can verify the correctness of our results. 

Our example, given the experience that we have 
accumulated along the past years, will refer to 
the O(N) a- model in two dimensions with lattice 
action 

1 



S 



2g L ^ 



(5) 



with <7 e S N - 1 and {dj) x = f x+fi - f x . The 
fields are related to their continuum version by 



a x = \/Z Z (t{x) . (6) 

All the renormalization factors have an ex- 
pansion of the form 



(=0 



n=Q 



•W In" 



(7) 



For example Zg and Z L are known up to four 
loops ||. 

For our purposes we made use of the non- 
perturbative determination of Zg given in || for 
the case N = 3. 

The Noether currents related to the O(N)- 
invariance, which is preserved also on the lattice, 
are 



and their singlet product is 



£# 6 (*)# ,6 (o) 

a . b 



(8) 



(9) 



In the continuum, using the OPE and averaging 
over a circle, we obtain 



|-Ji(*Wo(0) = 



9 



/d6 
— j x {r cos 0,r sin 6) • j (0, 



0) 



= [5 lO --9 o-] MS (0)W^°(r) + O(r 2 ) 

= gT^ s (Q)Wh°(r) + 0(r 2 ) (10) 

where T MS is the continuum energy-momentum 
tensor (EMT), which is exactly conserved, a prop- 
erty that is not shared by its lattice counterpart. 
It follows that 



x, t) dx = p, 



(11) 



where p is the momentum, so that, in a strip of 
size L, the expectation value between states of 
momentum p is 

(p\T^ s (0)\p) = p 
(p\p) L' 

The Wilson coefficient is given by 



(12) 



Wh\r) = 1 + 



5(N- 2) 



(13) 



TV -2 

47T 



.g( 7 + ln7r 2 A< 2 r 2 ) + 0(g 2 ). 



We have performed a Monte Carlo simulation on 
a lattice of size L x T — 128 x 256 with peri- 
odic boundary conditions at g^ 1 = 1.54, which 
corresponds to £ = 13.632(6). We measure: 
a) the correlation function in the one-particle sec- 
tor with momentum p 



C( P ,2t) = ±£\ 



S(p) 



-2tui(p). 



ip(x-y) 



- M»' ' (14) 

b) the correlation function of the out-of-diagonal 
entry of the lattice EMT with two one-particle 
operators 

Ti (p,2t) = 

= ~j- E((<9l<To,0 • 9o<To,o)((T-t,x ■ &t,y)) 



(15) 



xe ip{x-y) 

„ Z{p) (p\Tt o {0)\p) „ 
2w(p) (p|p) 

here we used the symmetric definition of the 
derivative 

Ix+ii f x — ii _ 



(16) 



c) the correlation between the products of two 
currents in the singlet sector and two one-particle 
operators 



1fiu(z;p, 2t) 



x,y 

Z{p) (p\3„,o-3 u Jp) 



ip(x-y) 



-2tui(p) 



2w(p) (p|p) 



(17) 
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■ p = 2jr/L 

■ p = 4ji/L 
p — 6jt/L 
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Figure 1. Angular average gj j l i(z;p,20)/2. 
The abscissa of coincident points have been split- 
ted. Notice that the rightmost points are not a 
faithful angular average because I has been mea- 
sured only on a square. 



Figure 2. Angular average gj j Xoi(z;p, 20)/2 di- 
vided by the theoretical prediction: g- (Wilson 
coefficient) -(energy-momentum tensor). Here in 
the Wilson function the first two leading logs have 
been summed up. 



Moreover we need a lattice version of the angular 
average 











(\z\-r-\ 






-i-w 











N(r) = £3(*) 



E(z) 



Let us now discuss the Monte Carlo data. 
In Fig. |l| wc plot the angular average 
g\lox{z;p, 20)/2 for various distances z and 
various momenta p — u2tt/L for n = 1,2,3. 



In Fig. H we plot the left-hand side of (10) 
divided by gTf$ s (0)Wq C '(r) renormalized at the 
scale ~p = 2.8519a. Here in the Wilson function 
the first two leading logs have been summed up. 
We find a large window of distances for which the 
OPE is verified, giving the correct matrix element 
independently of the external states. This win- 



dow seems to extend even at a distance of order 
£/2 which is presumably only a lucky accident. 

We also looked at the OPE directly on the 
lattice, i.e. using the Wilson coefficient com- 
puted in lattice perturbation theory resumming 
the first two leading logs, the expectation value 
Xio (p, 2t) , and the non-perturbative renormaliza- 
tion constant of the lattice EMT. In this case we 
do not observe a reasonable window where the 
OPE is verified. 

We wish to thank G. Martinelli and G.C. Rossi 
for useful discussions on this work. 
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